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We propose a variation of Bell inequalities for continuous variables that employs the Wigner 
function and Weyl symbols of operators in phase space. We present examples of Bell inequality 
violation which beat Cirel'son's bound. 
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Since its discovery in 1964 [1] the Bell inequality (BI) 
has triggered an enormous interest in the differences be- 
tween classical and quantum correlations. Bell inequal- 
ities are now commonly referred to as equations that 
relate correlation measurements which are fulfilled by 
any local hidden variable (LHV) theory, but are violated 
within the framework of quantum mechanics (QM). The 
original inequality was formulated for dichotomic vari- 
ables in spin systems. Clauscr, Home, Shimony, and Holt 
(CHSH) [2] presented a BI that was more amenable for 
experimental tests and is nowadays widely used. 

The original BI was inspired by the Einstein- Podolsky- 
Rosen paradox [3, 4] for continuous variables (CV). Bis 
for CV systems were first developed using dichotomic 
variables that have eigenvalues ±1 [5-7] . Recently, a 
new approach for CV Bell inequalities has been developed 
by Cavalcanti, Foster, Rcid, and Drummond (CFRD) 
[8, 9]. The CFRD inequality can be formulated for arbi- 
trary CV observables, but a multi-partite quantum sys- 
tem with five or more spatially separated sub-systems is 
usually required to obtain a violation. In this paper we 
offer a generalization of CFRD inequalities that employs 
the Wigner function W(q,p). The resulting Bis are con- 
ceptually different from those obtained using operator 
methods. Furthermore, the similarity of quantum phase 
space expectation values to those of classical statistical 
mechanics can provide new insight to the BI. 

The Wigner function is a quasi-probability distribution 
in phase space that is defined (in a ID system) by [10] 
W(q,p) = (2 7 rn)- 1 Smb[ / 5](g,p), with 



Smb[A](q,p)=2 / (q - q'\A\q + q') 



q')e 2tpq ' /h 



dq' (1) 



the Weyl symbol of an operator A and p the density 
matrix. W(q,p) appears in the expectation values of A, 



(A) 



dq dpW(q,p) Smb[A](q,p) 
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Phase space methods have been used to study specific 
implementations of the CHSH inequality for dichotomic 
variables [11-14]. Here we derive a BI that is based on 
Weyl symbols for a large class of observables A with di- 
chotomic symbols. As there is no general relation be- 



tween bounds of operators and bounds of their symbols, 
the resulting BI is conceptually different than other Bis. 

Local Hidden- Variable Theories — In the context of di- 
chotomic observables such as photon polarization, LHV 
theories are based on the following assumptions. 
LHV1: the possible measurement values of a LHV ob- 
servable B are given by the spectrum of the correspond- 
ing quantum observable B. 

LHV2: the expectation value of B is given by the expres- 
sion (-B)lhv = JdXf dbbp\(b), where the integration 
(or summation) variable b runs over the spectrum of B. 
Here, A represents all hidden variables in the LHV theory, 
and < p\(b) < 1 denotes the probability to measure the 
value b if the hidden variable takes the value A. Conser- 
vation of probability requires that J dX J dbp\(b) = 1. 
LHV3: For simultaneous measurements on two spatially 
separated systems, causality implies that the probabili- 
ties for the measurements on the two systems must be 
independent (i.e., they usually factorizc). 

In the context of CV observables, the assumption 
LHV1 is problematic. To see this, consider the case of a 
single ID particle. If we choose B to be the position 
q or momentum p of the particle, then the spectrum 
of B is the set of real numbers. Let us instead choose 
B = H = \{q 2 + p 2 ), which is the energy of a harmonic 
oscillator. If, in a LHV theory, q and p can take any 
real value, then the LHV spectrum of B should include 
all positive real numbers. However, the quantum me- 
chanical energy spectrum is discrete. Hence, if we tried 
to describe position and energy simultaneously, we could 
run into a contradiction about the spectra of the two 
observables. 

The origin of this contradiction is of course that, in 
quantum mechanics, H and q do not commute and thus 
cannot be measured simultaneously. On the other hand, 
in LHV models the observables must be commutative for 
at least a very large class of models [15]. The deriva- 
tion of our generalized BI is based on the assumption of 
commut at i vity. 

Re-derivation of the CHSH inequality. — Our method 
to derive Bis is related to the proof of the CHSH in- 
equality presented by CFRD [8] , which employs that the 
variance of a (generally complex) observable B must be 
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positive, |(-B)| 2 < (|i?| 2 ). The expectation values in this 
expression may either be evaluated within QM, then de- 
noted by (B), or in the framework of LHV models where 
we use the notation (-B)lhv- Within each theory this in- 
equality is always fulfilled. However, the maximum value 
of (|-B| 2 )lhv m all LHV theories provides an upper bound 
on local realism: if the quantum mechanical expectation 
value does not fulfill the inequality 

|(B)| 2 <max(|B| 2 ) L Hv, (3) 

Lrl V 

then the predictions of QM arc inconsistent with the as- 
sumptions behind LHV models. 

To derive the CHSH inequality we choose the observ- 
able B = X\X 2 + X\Y 2 + Y X X 2 -Y1Y2, with Xi, Yt four 
dichotomous observables (so that X? = Y 2 = 1) for 
two particles i = 1,2. In quantum mechanics one finds 
(B 2 ) = 4 + ([Yi,Xi] [X 2 ,Y 2 ]). CFRD then argue that 
in any LHV model the commutators must vanish, which 
leads to the CHSH inequality \{B)\ < 2. QM violates 
this inequality for a suitable choice of states and observ- 
ables. Because of the above-mentioned commutativity of 
many LHV models, this seems to be a suitable approach 
to deriving Bis. 

Bell inequalities in phase space. — In Ref. [8], CFRD 
used the method of commuting observables to find Bis 
for a more general form of the Bell observable B. Here 
we employ this approach to derive Bis for a system of 
n degrees of freedom (e.g., n one-dimensional particles) 
in phase space. We consider an operator B that is a 
function of position and momentum operators qi,fti (i = 
1, • • • n) and will establish an upper bound for (B) in LHV 
models by ignoring the commutators between position 
and momentum. We represent B by a sum of Weyl- 
ordered products, 

B= B ^-^ ' ' ' ^ 2 "- 1 « Weyl , (4) 

with complex coefficients B^ 1 ... Un . Wcyl ordering of a 
product of operators A\, A 2 , ■ ■ ■ corresponds to the com- 
pletely symmetric sum 

{MM ■ ■ • A fc )wcyi := Tj ^ Ap 2 ■ ■ ■ A Pk , (5) 

p 

where stands for the sum over all permutations of 
the k operators. Any operator that possesses a Taylor 
expansion can be brought into the form (4). For instance, 
the operator B = qp can also be written as B = | (qp + 
pq) + if, which corresponds to a single degree of freedom 
(n = 1) and B hl = 1, B ofi = i\ in Eq. (4). 

In phase space, the left-hand side (l.h.s.) of Eq. (3) can 
be evaluated using relation (2). The Weyl symbol of a 
Weyl-ordered operator B can be evaluated by replacing 
the operators q\,Pi by the respective phase space vari- 
ables, Smb[B](qi,pi) = B(qi,pi). For the operator (4), 



the symbol B(qi,pi) then has the explicit form 
B( qi , Pi )= B m-^Pl 2 ■■•C"- 1 < 2 "- (6) 

Mir" ■C2» 

To find the right-hand side (r.h.s.) of Eq. (3) we eval- 
uate (BB^) in QM and eliminate the contributions of 
all commutators between position and momentum oper- 
ators. In quantum phase space one has 




where the star product between two operator sym- 
bols captures the non-commutativity and non-locality of 
quantum observables and is given by [16-18] 

/*$(*) = / / d 2n yd 2n zf(y)g(z) (8) 

x exp — (jj-Jz + z-Jx + x- Jy^j . 

Here, x — (q,p) € R 2 " and J = [ j ], where I de- 
notes the n-dimcnsional identity matrix. 

We now employ the central assumption that the LHV 
upper bound can be found by eliminating all commuta- 
tors between position and momentum operators. In the 
Weyl symbol representation this can easily be accom- 
plished by taking the limit h — > in the star product 
B -k B* of Eq. (7). For smooth, h- independent /, g the 
star product has the expansion / * g(x) = f{x)g{x) + 
lH/> 9}( x ) + 0(^- 2 ); where the ?i-linear term is a Poisson 
bracket [16-18]. Consequently, B * B* in Eq. (7) is re- 
placed by the algebraic product \B\ 2 so that a tentative 
new BI is given by 

/n 
Y[d qj d Pj \B( qi , Pi )\ 2 W( qi , Pi ). (9) 

The r.h.s. of Eq. (9) does not yet provide the cor- 
rect upper bound for LHV theories because it depends 
on the Wigner function of QM. Intuitively one could fix 
this by replacing the quantum state by a suitable phase 
space distribution M / LHv(9ijPi) which is compatible with 
the assumptions of LHV models. However, general LHV 
models do not necessarily possess such a phase space dis- 
tribution. 

Instead, one has to consider suitable Bell operators for 
which the r.h.s. of Eq. (9) can be evaluated without re- 
ferring to a specific state. A simple but relevant example 
is the case when \B(qi,pi)\ 2 = \B \ 2 is constant. This 
happens if B(qi,pi) is a phase factor or a sign function, 
for instance. Because of the normalization of the Wigner 
function the BI then becomes 

KB^mi 2 <\B \ 2 . (10) 
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This is the main result of our paper. Wc remark that 
the dependence on the state has been removed here be- 
cause the integral of the Wigner function over the entire 
phase space equals to the probability to find any mea- 
surement result, which is unity both in QM and LHV 
theories. Consequently, |So| 2 can be interpreted as the 
upper bound for commutative LHV theories. 

We emphasize that the BI (10) is generally different 
from CFRD or CHSH type Bell inequalities. First, the 
upper bound is independent of the quantum state, while 
for the CFRD inequality the state dependence needs to 
be addressed. This can be a rather subtle issue: some en- 
tangled quantum states do have a positive Wigner func- 
tion [19] and positivity of the Wigner function is not suf- 
ficient to ensure consistency with LHV models [20]. In 
fact, Revzen et at [13] have shown that a dichotomic CV 
BI can be violated with a non-negative Wigner function. 
Second, the upjper bound is determined by the symbol of 
the operator B rather than its spectrum. This can be 
exploited to find new examples of BI violation. 

Example of Bell inequality violation. — We consider 
a system of two harmonic oscillators, which are prepared 
in the Bell state |^ Be u) - ^(|0) ® |1) - |1) ® |0)). 
We set H = 1 and measure lengths in units of the 
ground state width of the oscillators so that (q\m) = 

e~ V H m {q) I (\/2 m ml y/n) , with H m (q) the Hermite poly- 
nomials. If wc introduce a phase space variable x := 
(q, p) with x 2 = q 2 + p 2 , and employ center-of-mass vari- 



ables (5x := (xi — X2)/-\/2 and x c := (xi 
Wigner function of I'i/'Beii) is given by 



X2 



Wbo1i(x c ,<5x) = 7T 



^ (2<5x 2 - 1) 



)/y / 2, the 



(11) 



This corresponds to the product Woo{xc)Wu(6x.) of the 
ground state in the center-of-mass coordinates and first 
excited state in the relative coordinates. To exploit the 
negativity of the Wigner function we consider a Bell op- 
erator that has the dichotomic Weyl symbol £?(xi.X2) = 
sgn(2<5x 2 — l), which implies that the bound |-Bo| 2 m 
BI (10) is unity. We remark that the assumption that 
the Weyl symbol is dichotomic is different from the as- 
sumption that the spectrum of an operator is dichotomic. 
There is no general relationship between the bounds for 
an operator and the bounds for its symbol. 

The l.h.s. of BI (10) can be evaluated using Eq. (2), 
so that the BI turns into 



\(B) 



d 2 6x <i 2 x c W(6x, x c ) sgn (2 5x 2 - l) 



< 1 



(12) 

For the Bell state (11) we find \(B)\ = 4» -1 « 1-426, 
so that the generalized BI is violated in quantum me- 
chanics. This violation is slightly larger than Cirel'son's 
bound v* f° r the ratio between the maximal quantum 
mechanical value of \(B)\ and the bound of the CHSH 
inequality for LHV theories [21]. 



Generalization. - The previous example of BI vi- 
olation exhibits two special features: (i) in center-of- 
mass variables, the Wigner function W takes a product 
form. The negativity of W depends on just the relative 
phase space variable <Sx. (ii) The Bell operator's sym- 
bol B(Sx) exactly cancels the negativity of W, so that 
(B) = J d 2 Sx\W(5x)\. These features can be readily 
generalized. 

Consider a Wigner function W(x), with x G R 2n , that 
is negative on a measurable subset £_ C R 2 ™ and posi- 
tive on its complement £+. We introduce characteristic 
symbols x±(x) that are unity for x £ £± and zero oth- 
erwise. The symbol of the Bell operator is then given by 
B(x) = x+( x ) — X-( x )- The general relation between an 
operator on Hilbcrt space, % = L 2 (K n , C), and its Weyl 
symbol is given by 

A = y"d 2 "xA(x)SmbL4](x), (13) 

with the quantizer [18] defined via its Dirac kernel 
(q'\A( X )\q") = (nK)- n 5 (q' + q" - 2q) e *f(«'-«"). (14) 

The two operators corresponding to the characteristic 
symbols are therefore given by \± = J £ d 2n x A(x). To- 
gether, they form a partition of unity, 1 = x+ + X-^ ancl 
commute. The non-local character of the star product 
(8) causes \+ * X- an d as a result the Bell operator 
fulfills B 2 ^ t. This implies that B is not a dichotomic 
operator, even though its symbol is dichotomic, in the 
sense that B(x) 2 = 1. 

As a concrete example, we consider the special case 
that £- consists of a disk of (dimensionless) radius R 
around the origin of a 2D phase space. Since X- ( x ) G 
i 2 (R 2 ,C) the corresponding operator x- 1S Hilbert- 
Schmidt and so has a discrete spectrum. Furthermore, 
because of its spherical symmetry commutes with the 
Hamiltonian of the harmonic oscillator, so that the eigen- 
states of X- ar e the harmonic energy eigenstates |m). 

The eigenvalues X m (R) of X- are conveniently com- 
puted from the fact that they are the expectation values 
with respect to the state \m) [22] 

\ m {R) = ^ / d 2x X-« Smb [\m)(m\] (x) . (15) 

The integral above is the evaluation of the trace 
Trx_|m)(m| via the corresponding Weyl symbol. Us- 
ing Smb[|m)(m|] (x) = 2(-l) m i,„(2x 2 ) exp(2x 2 ) where 
L m is the Laguerre function gives 



X m (R)= / (-l) m L m (r 2 )e- r dr 2 . (16) 
Jo 

By replacing the Laguerre polynomials with their gener- 
ating function (Eq. (22.9.15) of Ref. [23]), we can derive 
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1 - 2 A, 




FIG. 1: BI bound |1 — 2A m | < 1 for a system prepared in state 
\m) and dichotomic operator symbols corresponding to a disk 
of radius R in phase space. Shown are the cases R = l/v2 
(circles), R = 3.5 (squares), and R = 5.5 (diamonds). 



In conclusion, we have proposed generalized Bis for 
which the bound is determined by the Weyl symbol of 
the Bell operator. Examples of states for which the BI 
is strongly violated have been presented for bi-partite 
systems. Extensions of this work for multi-partite or 
interacting systems may reveal further insight into the 
differences between classical and quantum mechanics. 
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FIG. 2: BI violation of \B\ = / d 2 x|W m (x)| < 1, with W„ 
the Wigner function of state |m). 



a generating function for the eigenvalues of x_, 

GWHt-ir 1 ^-!)- (17) 



The mth eigenvalue A m is given by the mth coefficient of 
the Taylor expansion of G(t) around t = 0. 

For the Bell operator 1 — 2£- and the relative coordi- 
nate prepared in an eigenstate of \- > the BI then turns 
into the condition |1 — 2A m | 2 < 1. In Fig. 1 we display 
the eigenvalues 1 — 2 A function of the excitation 

number m for different choices of R. The example of 
the Bell state above corresponds to the case m = 1 and 
R = 1 1 v2- For larger R and higher excitation numbers, 
a similar degree of BI violation can be obtained. 

To increase the degree of BI violation, one can choose 
£- to agree with the area in phase space where the 
Wigner function W m of state |m) is negative, so that 
(B) = J cPx |W m (x)|. The result for this expression for 
m < 30 is shown in Fig. 2. BI violations that are much 
larger than in the CHSH case are possible. However, 
because of the normalization and exponential decay of 
the Wigner function we conjecture that the maximal BI 
violation for arbitrary m will be bounded. 
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